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We analyze dynamics of 3D coreless vortices in superfiuid films covering porous substrates. The 
3D vortex dynamics is derived from the 2D dynamics of the film. The motion of a 3D vortex is 
a sequence of jumps between neighboring substrate cells, which can be described, nevertheless, in 
terms of quasi-continuous motion with average vortex velocity. The vortex velocity is derived from 
the dissociation rate of vortex-antivortex pairs in a 2D film, which was developed in the past on the 
basis of the Kosterlitz-Thouless theory. The theory explains the rotation-induced dissipation peak 
in torsion-oscillator experiments on 4 He films on rotating porous substrates and can be used in the 
analysis of other phenomena related to vortex motion in films on porous substrates. 

PACS-numbers: 67.40.Vs, 67.57.Fg, 67.70.+n, 67.40.Rp 

I. INTRODUCTION 

Superfiuid 4 He films adsorbed in porous media is an actual topic in physics of superfluidity Studying of this 
system gives a unique possibility to investigate the interplay between 2D and 3D physics, especially the character of 
the transition to the superfiuid state. On one side, torsion-oscillator experiments reveal the dissipation peak near the 
temperature of the superfiuid onset T c , which is predicted by the dynamical theory of vortex-antivortex pairs [2J based 
on the Kosterlitz-Thouless theory for 2D films On another side, they found that in films on porous substrates 
the superfiuid density critical index ~ 2/3 [l[ of the 3D system and the sharp cusp of the specific heat Q at T c are 
similar to those near the A transition of the bulk 4 He. 

An important insight into physics of superfiuid films in porous media is provided by torsion-oscillator experiments 
with rotating substrates. The porous substrate is usually modeled with a "jungle gym" structure [!, [H, H HI : a 3D 
cubic lattice of intersecting cylinders of diameter a with period I (Fig[JJ). Multiple connectivity of superfiuid films 
in porous media allows a variety of vortex configurations, and probably most important from them is a coreless, or 
pore 3D vortex, which is just a flow around the vortex pores having nonzero circulation. Due to the presence of a 
new type of topological defects, one could expect an essential difference in the response between the plane film and 
the porous-medium film under rotation. This expectation was confirmed by torsional-oscillator experiments, which 
revealed a rotation-induced peak in dissipation (inverse quality factor) as a function of temperature. The additional 
peak was shifted from the stationary (static) peak that was observed without rotation [l(| . Double-peak structure 
essentially differs from the case of the plane film, where the only effect of rotation was to broaden the stationary peak 
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A semi-empirical interpretation of the rotation-induced dissipation peak was suggested in Ref. [LQ. It is clear that 
the rotation can affect dissipation via rotation contribution to the velocity field. So it is a nonlinear correction to the 
response. Instead of the derivation of such a correction from the theory, the authors of Ref. [1(3 used the data on the 
nonlinear response taken from the independent experiment on large-amplitude torsion-oscillations. This provided a 
qualitative explanation of the rotation peak, and even of some quantitative features of it, but could not pretend to 
be a full theory of the effect. The present work suggests really a theory of the rotation peak deriving the parameters 
of the peak from the parameters of the film and the substrate. The key role in our scenario is played by the 3D 
coreless vortices. We derived dynamics of these vortices from dynamics of 2D films adsorbed in porous media. It was 
suggested [l(| that motion of the coreless 3D vortices occurs in a creeping manner by jumping from cell to cell. These 
jumps are related to dissociation of the vortex-antivortex pair on one side of a rod separating different pores, with 
subsequent annihilation of the pair on another side of a rod. The result of the jump is a shift of velocity circulation 
to a neighboring pore. Though our analysis was focused on application to the torsion-oscillator experiments, it is 
valid for description 3D vortex motion in many other cases, at least in those, where the vortex moves to distances 
much larger than the average period of the substrate structure. An important example is steady vortex motion in 
zero-frequency experiments. 

In Sec. [IT] we describe how dynamics of 3D coreless vortices is connected with dynamics of 2D films covering the 
substrate. This provides a bridge between the effective-continuous-medium 3D description and dynamics of 2D films. 
Section [HI] reviews the theory of dissociation of vortex-antivortex pairs on the basis of the Kosterlitz-Thouless theory. 
In Sec. HVI the pair-dissociation rate is analyzed near the critical temperature, where there is an analytical solution of 
Kosterlitz's recursion equations. The theory is applied to the analysis of the torsion-oscillation experiments in Sec.fVl 



FIG. 1: (Color online) Vortex creep in the jungle gym structure: the coreless vortex (dashed line) crosses a cylinder between 
two cells via dissociation of the vortex-antivortex pair. The arrowed curved lines show circulation around the cells where the 
vortex line is located. 

The rotation-induced dissipation peak directly follows from the theory, and its position, shape and even height are in 
a good agreement with the experiment. Section IVT1 is devoted to conclusions. 

II. CORELESS VORTICES: FROM 2D TO 3D VORTEX DYNAMICS 

In a continuous medium a vortex is a topological defect with nonzero circulation around the vortex axis. There 
is an area around the vortex line with a suppressed order parameter, which is called vortex core. But topology of 
a porous medium allows circulation of superfiuid velocity around a pore without suppression of the order parameter 
anywhere inside the superfiuid film. This leads to the concept of a coreless vortex. The vortex "line" in this case is 
not a line at all; this is a chain of the jungle gym structure cells with nonzero circulation around them. Schematically 
the coreless vortex is depicted in Fig. Q] This coreless structure of a vortex rules out usual type of vortex motion in 
a continuous medium simply because the coreless vortex has no continuous coordinate: its position is discrete and is 
determined by a cell with nonzero circulation around it. The only way for the vortex to move is to jump from cell to 
cell. They call such a type of vortex motion vortex creep. But any jump is in fact a process in time (whatever short) 
during which the vortex line inevitably cross a rod covered by a superfiuid film. So during the jump the "coreless" 
vortex does have cores: at the place where it enters the rod and at the place where it comes out from the rod (Fig.QJ. 
These two cores together form a 2D vortex-antivortex pair, which should grow, dissociate, and eventually annihilate 
on the other side of the rod. This process lead to a discrete shift of the vortex line to a neighboring cell. For better 
illustration a two-dimensional picture of this process is shown in Fig. [2j It shows two cells: upper (u) and lower (1). 
Before the process (Fig. (5^) there is circulation k = h/m,He around the lower cell as shown in the figure. Provided 
that there is no other coreless vortices nearby, the same circulation exists around any path inclusive of the cell (1) as 
shown for the path around the both two cells. The transient process of the "jump" from the cell (1) to the cell (u) is 
shown in Fig. [2)3. The vortex-antivortex pair is present in the film between two pores. This makes circulation around 
any of the two cells undefined: it depends on whether the path goes outside or inside the pair. Only circulation around 
the two cells together remains equal to k. Figure shows the state after the process: the coreless vortex is now 
located in cell (u). Though this scenario is shown in the plane picture, it is directly applicable to the 3D jungle gym 
structure with the 2D vortex and antivortex moving around a cylindric rod. In reality creation and dissociation of 
vortex-antivortex pairs is a stochastic process, which is determined by the average dissociation rate R(T, v s ) (number 
of dissociation events per second and per unit area of a film) , which depends on temperature T and, most important 
for us, on average superfiuid velocity v s in the film. In our case the superfiuid velocity field consists of two parts: 
v s = v c + V s , where v c originates from a circular flow around the vortex line at rest and V s is a transport superfiuid 
velocity with respect to a moving substrate. In Fig. [2] the velocity V s is directed along the axis x. So it it is added 
to the circular velocity v c above the cell with circulation k and is subtracted from v c below this cell. Therefore the 
jumps of the vortex line up and down are unbalanced and result in some average drift (creep) of the vortex line up the 
picture in Fig. [1] This jump-like process can, nevertheless, be described with an average vortex velocity determined 
as 

f) /? 

V L = Al(R(v c + V S ,T) - R(v c - V S ,T)) « 2AI — V s , (1) 

ov s 

where I is the period of the jungle gym structure, and A is the area of the rod separating neighboring pores. According 
to this scenario of vortex motion the vortex velocity is strictly normal to the velocity V s . 

Let us compare the relation ([1]) with the general relation connecting the vortex velocity of the vortex line with the 



3 




VAP 





=K 



a) 



b) 



c) 



FIG. 2: Vortex creep. The superfluid moves with velocity V a along the axis x, and the vortex creeps along the axis y. (a) The 
state before the jump. There is circulation k around the lower cell (1) as it is shown in picture. The same circulation exists 
around any path inclusive cell (1). The circulation around two cells (1) and (u) is shown. The coreless vortex is located at the 
cell (1). (b) The vortex on the way between cells (1) and (u). There is the vortex-antivortex pair (VAP) in the film. Circulations 
around the paths inclusive only of cell (1) or (u) are not defined, while for the path around the two cells circulation is equal 
to k as before, (c) The state after the jump. The coreless vortex is located at the cell (u). No circulation around the cell (1) 
anymore, but there is circulation around the cell (u), or around the two of them. 



normal and the superfluid velocities [1 31 ] : 

V L = V s + az x (V„ - V s ) - a'(V„ - V s ) . 



(2) 



In the problem under consideration we can assume that the normal component is clamped to the porous-glass substrate 
oscillating with the velocity V g . So their velocities coincide: V„ = V g . In fact equation |T]) is written for the system 
moving with the substrate, where V„ = V s = 0. Eventually Eqs. |T]) and @ agree if 
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Note that the condition 1 + a' — providing that vortices move normally to the superfluid motion means the absence 
of the "effective" Magnus force [14[, which is defined as the term oc [i X Vj;] is the balance of forces on the vortex. In 
superconductors this leads to the total absence of the Hall effect. Absence of the effective Magnus force is typical for 
lattice systems [lBj . in contrast to uniform continuous media with Galilean invariance. The crossover between these 
two cases was recently studied by numerically solving the Gross-Pitaevskii equation flfjj ]. 

So we have demonstrated that the creep of coreless vortices, which is realized via sequences of discrete jumps from 
cell to cell, can be described in the terms usually used for 3D vortices moving in a continuous medium. Still the 
parameters of this 3D "effective-medium" description must be determined within the theory of 2D films covering the 
multi-connected substrate. The crucial parameter to be determined is the derivative dR/dv s of the dissociation rate R. 



III. RATE OF PAIR DISSOCIATION 

The thermally activated dissociation of vortex-antivortex pair was analyzed by Ambegaokar et al. [2| (see also 
references to later works in Ref. ilTi) on the basis of the Kosterlitz-Thouless theory. They considered superfluid films 
on plane substrates, while in our case films cover cylindrical surfaces. But as we shall see below, the relevant scale 
(size of the pair at the saddle point) is small compared to the substrate curvature and the curvature may be ignored. 

The pair dissociation is accompanied by overcoming the potential barrier. The barrier corresponds to the saddle 
point of the vortex-pair energy as a function of the radius-vector r connecting vortex with antivortex, 

TTt \ P s0K2 f T d2lC I *\ fA\ 

u (r, v s ) = — — / ~, r - p s0 Kr ■ (v s x z) . (4) 
2tt Jao £{r)r 
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Here p s o is the bare superfluid density, e is Kosterlitz's static scale-dependent dielectric constant determined from the 
integral equation 
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ro is the core radius of the 2D vortex, yo = er E °l T , and Eq is the energy of the vortex core. The dielectric constant 
e(r) takes into account screening of interaction between a vortex and an antivortex at distance r by pairs of smaller 
size. In the limit r — > oo the dielectric constant determines the ratio of the bare and the renormalized superfluid 
densities: p s — p s o/e(oo). 

The integral equation (O can be reduced to Kosterlitz's recursion equations: 
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Here / = ln(r/ro), the dielectric constant is replaced by the ratio e(7) = K(0)/K(l), where K(0) = p S QK 2 / '4tv 2 ksT is 
the bare Kosterlitz -Thouless coupling constant related to the bare superfluid density p s0 , and 
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is the rescaled activity. 

Solving the Fokker-Planck equation for the distribution of vortex-antivortex pairs, Ambegaokar et al. Q obtained 
the following expression for the dissociation rate [see Eq. (4.6) in their article] 
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R=^y 2 (l s )exp(2nK(l s )). 



(8) 



Here D is the vortex diffusion coefficient, r s is the saddle-point value of the pair size r, and l s = r s /ro is its logarithm. 
The saddle point for the effective potential (H]) is reached when r is perpendicular to v s and r s satisfies the condition 
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If the velocity v s is small the values of l s and r s are large, and this condition yields that 
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In order to find the mutual friction parameter a from Eq. ([3]), we need the derivative of the dissociation rate with 
respect to the velocity v s . Using Eqs. (O and © we obtain: 
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In summary, Eqs |TJ), |3|) and show how the creep motion of the coreless vortices in a porous medium is 
described in terms of parameters determining motion of 3D quantum vortices in continuous media. This description 
can be used for various problems related to vortex motion in porous media. 



IV. PAIR DISSOCIATION NEAR THE CRITICAL POINT 



For understanding the nature of the rotation dissipation peak we need to study the temperature dependence of 
the dissociation rate at temperatures close to the critical one. At these temperatures Kosterlitz's recursion equations 
have an analytical solution 0]. One can introduce a small x(l) = n[K(l) — K c (oo)}, where K c (l) yields values of K(l) 
at the critical point and at I — > 00 K c (oo) = 2/tt. Then the recursion relations can be written as 
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Their solution is 

m / , , -i %o \ x cosh(x 00 /) + x 00 siah(x 00 l) 

x{l) = coth [Xool + coth =2:00 —7 7T- T7 tt, (13) 

V ^oo / a;osinh(ar 00 <) + x 00 cosh(a; 00 <J 

4^(0 = a^csch far^i + coth" 1 = . - — — . (14) 

Here xo = x(0) and Xoo = x(oo). The solution satisfies the condition 

xl=x(l) 2 -[^y(l)] 2 = x 2 - (47n/ ) 2 . (15) 
At the critical point Xoo — and Eqs. (fl~3|) and (I14[) become 
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According to Eq. (|15[) . at the critical point the parameters xq and yo satisfy the relation 
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As usually assumed in the 2D vortex dynamics, the bare superfiuid density does not vary near the critical point. 
Then, since K(0) oc 1/T and yo — e~ E °/ T , their dependence on the relative temperature t = (T c — T)/T c is 

K(0) « K C (0)(1 +t), y « yoc (l - ^tj . (17) 

Then Eq. yields that at t > x^ = 2b^/i with b = 167re~ B °/ Tc [1 + 2tt(1 + So/T^e-^"/ 7 ^]. This leads to the 
square-root cusp 

Ps {T) = Ps (T c )(l + bVt) (18) 

in the critical behavior of the renormalized superfiuid density, which was revealed by Nelson and Kosterlitz (l8j with 
numerical calculations. 

Using all these relations together with the assumption that v s so low that at the saddle point l s 3> xoc, we obtain 
linear dependence of the dissociation rate on t: 

^ = »v_ 7() , (19) 

avs 2 Kr s<-s 

where 

46 2 

7 = 2b% + — - 1 . 

V. TORSIONAL OSCILLATIONS OF ROTATING POROUS SUBSTRATE AND COMPARISON WITH 

THE OBSERVED ROTATION DISSIPATION PEAK 

Let us now apply the theory to oscillatory motion of the substrate superimposed on its steady rotation with the 
angular velocity SI. So the substrate velocity field is [SI x r] + V s . Only the oscillatory component V g oc e _iwt is 
important for us. Using the Euler equation 

^ + [20x V L ] = -V/i (20) 

and Eq. ^ with a' = — 1 one obtains for an oscillatory components of the velocities (the chemical potential [i is not 
relevant for the azimuthal motion): 



(21) 
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This relation determines the drag of the superfluid component by the oscillating substrate. Analyzing now the balance 
of forces for the torsional resonator, as was done many times in the past, one obtains the following contribution to 
the inverse quality factor of the torsion oscillator for the slow rotation Via <gC uj: 

where M is the total mass of the torsional oscillator, V is the total 3D volume (including pores), and p s s is the 
effective 3D superfluid mass density in the porous-glass substrate of the volume V, which is connected with the 2D 
superfluid density p s of the film by the relation 

Ps* = y ■ (23) 

Here A to t is the total area of the film. For the "jungle gym" structure with cell size I and the rod diameter a 

^£ « — . (24) 

On the other hand, in the theory of 2D superfluid films starting from Ref. [H they describe the drag of the superfluid 
component introducing the "dynamical dielectric constant" e(uj) (one should not mix it up with Kosterlitz's static 
dielectric constant e introduced in Sec. : 



Then 
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with the imaginary part of the inverse dielectric constant equal to 

1 2Vla 4AIQ dR , s 

Im- = 5T- • (27) 

e uj uj ovs 

Substituting Eqs. (fig)) and ([!§)) into Eqs. and (|2"7|) . we see that the square- root cusp in the superfluid density 
is crucial for its temperature dependence in the critical area and for existence of the rotation dissipation peak: 

A Q-i = |fl = AQ-\T C )(1 + bV~t- <yt) . (28) 

uj M ov s 

The factor 7 oc l c is expected to be large, but the the square-root cusp is more essential at small t, and the inverse 
quality factor has a maximum at t = b 2 /A^ 2 rather close to the critical point. 

The linear approximation used for derivation of Eqs. (fT3"|) and (fl"4"|) is more or less truthful only for \T — T c \ 
not exceeding 0.005 K. Though the dissipation maximum at T w 0.6234 K is in this interval (see below), the low 
temperature side of the dissipation peak is not, and the peak width cannot be determined using the analytical formulas 
(H~5]) and (Tj"4"| . Thus though the analytical expression (|2"8|) qualitatively explains the observed dissipation peak itself, it 
is not sufficient for its quantitative description. A more accurate quantitative analysis required numerical calculations. 

Let us now gather from Refs. Q and [l(| all quantitative data needed for comparison with the theory. The porous- 
glass substrate can be modeled with the jungle-gym structure with the diameter of rods a sa 1 /im and the structure 
period / « 2.5 /im. Then the circulation velocity around the pore is estimated as v s = n/Al « 1 cm/sec. According 
to Ref. [l(J, the transition temperature is T c w 0.628 K. The areal density of superfluid component at the critical 
temperature (jump of density) can be evaluated from the Kosterlitz-Thouless relation 

87rfcsT c — q 2 

p s = 5 — = 2. 179 1 x 10 a g/cm . 

K 

We also need the bare areal superfluid density p sa . It can be calculated from dependence of the transition temperature 
on the thickness: According to Ref. 0, the coverage of the substrate is about 33 x 10~ 10 mol/cm 2 . Then p S Q w 
1.32 x 10~ 8 g/cm 2 , and the bare Kosterlitz -Thouless coupling constant K(Q) is 
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FIG. 3: Comparison of the experimental and the theoretical inverse quality factor of torsional oscillations. The right peak is 
the static peak, which is present without rotation. The group of the left peaks is rotation-induced ones obtained experimentally 
for various £1 and scaled (divided) by fl. The fact that all these curves collapse well on the same curve proves linearity on 
fi. The dashed line shows the theoretical inverse quality factor. The vertical scale of it was fitted so that the height of the 
theoretical peak coincided with the experimental rotation peak, but the width and the position of the theoretical peak were 
calculated without any fitting. The theoretical curve ends with the critical point. The experimental superfluid density scaled 
by its zero temperature value is is also displayed (the left vertical axis). 

Quantity b entering relations (|13I) - ([T4")) can be obtained from the width ~ 0.02 K of the experimental curve describing 
static dissipation peak of paper [lOj]. It yields b about 3.6975. Furthermore, the diffusion coefficient is D = 2.8 x 10~ 7 
cm 2 /s, and the core radius of vortices on film can be estimated as ao ~ 25 x 10 -8 cm. 

We can now find the pair size r s at the saddle point from Eq. (fT0|) : r s w 2.98 x 10 -5 cm. This corresponds 
to l s — hx(r s /ro) — 4.78. On one hand, these values are large enough in order to justify our assumption of large 
l s compared to xq c . On the other hand, we see that r s << a, which justifies our consideration of vortex-antivortex 
unbinding using the theory for plane films. 

We also need the value of A tot p s (T)/M in Eq. ([2fJL This is the ratio of the temperature dependent superfluid 
mass to the mass of the empty cell. According to Ref. |9| it is approximately equal to 10 -5 . For determination of the 
temperature dependent superfluid mass p s (T) we have built the extrapolating function using the experimental results 
obtained from the measuring of the shift of the oscillations period (see Fig. 2 from Ref. KLOh . 

In order to calculate K(l, T) and y(l, T) numerically we use procedure proposed in Ref. Qjl (see Appendix A there). 
In the temperature interval from T c to the temperature T « 0.6234 K corresponding to the dissipation maximum 
the analytical expressions based on the linear approximation are more or less truthful. But in the low temperature 
region T < 0.6234 K we performed the numerical calculation using the Mathematica program. Following Ref. [IH, 
for any temperature \T — T c \ and corresponding x we choose a value Iq so that x^Iq = tt/2. Under this choice both 
y(lo,T) and deviation of K(Iq,T) from 2/ir are rather small and the analytical expressions |T5j) and (fT4l are still 
good enough for evaluation of K(l,T) and y(l,T). Further we take K(Iq,T) and y(lo,T) as initial conditions for 
numerical integration (with respect to variable /) of the recursive Kosterlitz-Thouless relations ([6]). In this way we are 
able to restore K(l s ,T) and y(l s ,T) at saddle point l s and further to find dR/dv s and eventually AQ -1 in the whole 
temperature interval T < T c . In Fig.?? we plot the theoretical AQ _1 (T) (dashed line) together with the experimental 
data of Fukuda et al. [l(| ■ The observed rotation induced dissipation peak (the left peak) was scaled by the angular 
velocity reducing it to the value measured at Q = 6.28 rad/s. The plots for different values of ft collapse on the 
same curve (see Fig. 2 in Ref. 03), which proves linear dependence of rotation induced dissipation on f2. For better 
comparison of the peak shape we fit the theoretical height of the peak to the experimental one. But in fact the values 
of AQ _1 (T) at the maxima do not differ essentially: At the angular velocity ft = 6.28 rad/s they are 3.54 * 10~ 8 
in the theory and 2.4 * 10~ 8 in the experiment (i.e., about 70 % from the theoretical value). Keeping in mind that 
nice agreement in the peak shape and position was achieved without any additional fitting the agreement is really 
satisfactory. We conclude the quantitative analysis of the torsional-oscillation experiment with the estimation of the 
creep velocity Vl of the pore vortex using Eq. (JTJ. For V s 1 cm/sec, which is typical in the torsional oscillation 
experiments, we get Vl ~ 0.7 cm/sec. For the frequency v = 477 Hz in the experiment, this corresponds to vortex 
displacements about 1.5 x 10 -3 cm, which is about 6 structure periods. This looks rather satisfactory for our scenario 
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reducing sequences of jumps between neighboring cells to quasi-continuous vortex motion. 



VI. CONCLUSION 



We suggested the theory of motion of 3D coreless vortices through a 2D film covering a porous substrate. Dynamics 
of 3D vortices is derived from dynamics of vortex-antivortex pairs in the 2D film. The 3D vortices move jumping 
from cell to cell of the substrate structure, but it can be described in terms of the average vortex velocity like vortex 
motion in continuous media. Frequency of jumps, and vortex velocity correspondingly, is determined by dissociation 
rate of vortex-antivortex pairs, which is known from dynamics of 2D superfluid films based on the Kosterlitz-Thouless 
theory. We calculated dissipation intensity and its temperature dependence analytically and numerically. The theory 
is compared with the experiments on torsional oscillations of 4 He superfluid film adsorbed on a rotating porous-glass 
substrate. We explain the second (rotation-induced) dissipation peak on the temperature dependence, which was 
revealed in these experiments. Quantitative comparison between theory and experiment looks satisfactory, especially 
for the shape and the position of the peak. 

Though we focused on application of our theory to torsional-oscillation experiments in He films on porous-glass 
substrates, we believe that the theory has a much wider area of possible applications. The concept of quasi-continuous 
vortex motion in porous media with parameters determined from dynamics of 2D superfluid films should be applicable 
not only to straight vortices induced by steady rotation. For example, one may apply the theory also to vortex rings 
and to the vortex tangle if they can be created in porous media. The theory can also be used for the analysis of 
steady vortex motion in various situations, e.g. in the process of heat transfer. 



ACKNOWLEDGMENTS 



Authors are grateful to M. Kubota for numerous discussions and providing necessary experimental data. S. K. 
N. thanks the Racah Institute of Physics of the Hebrew University of Jerusalem for hospitality and support and 
acknowledges a partial support by the grants N 05-08-01375 and 07-02-01124 of RFBR and the grant NSH-6749.2006.8 
of the state support of leading scientific schools by the President of the Russian Federation. . 



1] J. D. Reppy, J. Low Temp. Phys. 87 205 (1992) and references therein. 

2] V. Ambegaokar, B.I. Halperin, D.R. Nelson, and E.D. Siggia, Phys. Rev. B 21, 1806 (1980). 
3] J. M. Kosterlitz and D. J. Thouless, J. Phys. C6, 1181 (1973). 
4] S.Q. Murphy and J.D. Reppy, Physica B 165, Part 1, 547 (1990). 

5] J. Machta and R. A. Guyer, Phys. Rev. Lett. 60, 2054 (1988); J. Low Temp. Phys. 74, 231 (1989). 
6] T. Minoguchi and Y. Nagaoka, Progr. Theor. Phys. 80, 397 (1988). 
7] G. A. Williams, Phys. Rev. Lett. 82, 1201 (1999); ibid. 68, 2054 (1992). 
3] T. Obata and M. Kubota, Phys. Rev. B 66,140506 (2002). 

9] M. Fukuda, K. Ooyama, T. Obata, V. Kovacik, and M. Kubota, J Low Temp. Phys. 113, 417 (1998); ibid. 113, 423 
(1998). 

d] M. Fukuda, M.K. Zalalutdinov, V. Kovacik, T. Minoguchi, T. Obata, M. Kubota, and E.B. Sonin, Phys Rev B 71, 212502 
(2005). 

1] P. W. Adams and W. I. Glaberson, Phys. Rev. B 35, 4633 (1987). 

2] D. J Bishop and J. D Reppy, Phys. Rev. Lett. 40, 1727 (1978); Phys. Rev. B 22, 5171 (1980). 
3] E.B. Sonin, Rev. Mod. Phys. 59, 87 (1987). 

4] E.B. Sonin, in: Vortices in Unconventional Superconductors and SuperHuids, edited by R.P. Huebener, N. Schopohl, and 
G.E. Volovik, (Springer-Verlag, Berlin, 2002), pp. 119-145. 
[15] E.B. Sonin, Phys. Rev. B 55, 485 (1997), sec. V. 
[16] Z. Geese and S. Khlebnikov, Phys. Rev. B 72, 054525 (2005). 

[17] K. A. Gillis, S. Voltz, and J. M. Mochel, J. Low Temp.Phys. 61, 172 (1985); S. Giorgini and R. M. Bowley, ibid. 102, 171, 

(1996); R. M. Bowley, A. Kirk, and P. J. King, ibid. 88, 73 (1992). 
[18] D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 1201 (1977). 



